Surface State Stark Shift in a Scanning Tunneling Microscope 
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We report a quantitative low-temperature Scanning Tunneling Spectroscopy (STS) study on the 
Ag(lll) surface state over an unprecedented range of currents (50 pA to 6 /iA) through which we 
can tune the electric field in the tunnel junction of the microscope. We show that in STS a sizeable 
Stark effect causes a shift of the surface state binding energy Eq. Data taken are reproduced by 
a one-dimensional potential model calculation, and are found to yield a Stark-free energy Ea in 
agreement with recent state-of-the-art Photoemission Spectroscopy measurements. 

PACS numbers: 73.20.-r,68.37.Ef 



Surface states have been investigated intensely for over 
two decades of surface science. These states, which are 
trapped between the surface barrier potential and a band 
gap in the crystal, are an experimental realization of a 
quasi-two-dimensional electron gas with a characteristic 
dispersion of E(k) — Eq oc fc 2 , where Eq is the lower 
edge of the energy band. Scanning Tunneling Microscopy 
(STM) and Spectroscopy (STS) have provided a new way 
to study surface states through local measurements per- 
formed on the atomic scale. The first imaging of Shockley 
surface states on the (111) facet of the noble metals Ag, 
Cu and Au 1], has motivated a great deal of STM and 
STS studies on these systems. Since then, many aspects 
have been elucidated such as surface-state-mediated in- 
teractions 0, and lateral confinement to nano-cavities 
0. The lifetime r of these states at Eq, as determined 
by STS 0, IE 01 also agree remarkably well with re- 
cent state-of-the-art angle-resolved Photoemission Spec- 
troscopy (PES) measurements 0. Despite this success, 
a worrying discrepancy between STS and PES still sub- 
sists in noble metals concerning Eq, as -Eo(STS) is found 
to lie lower than i?o(PES) typically by —5 to —20 meV 
110. 

This last point is still an open question. A possible 
explanation would be that unlike STS, where the local 
character of the measurement allows to access a defect- 
free region of the surface, a PES measurement is defect- 
sensitive as it integrates over a large sample area of about 
1 mm 2 . And defects should, in principle, shift i? (PES) 
upwards in energy . However defects also broaden the 
PES line width (r(PES)= H/t), and the recent agree- 
ment found for r between STS and PES appears to in- 
validate this scenario. 

In this Letter, we present STS measurements on the 
surface state of Ag(lll) which elucidate this discrepancy. 
The presence of an electric field between the STM tip and 
the sample surface is known to produce a Stark shift in 
the STS spectra of field emission resonance states , or 
a more dramatic effect such as band bending in semicon- 
ductors With this in mind, we have tracked down 
the influence of the electric field on the STS spectrum of 



the Ag(lll) surface state. The tunneling current J in a 
STM junction can be written as 

I oc exp(-1.025V^ d) (1) 

where </> is the apparent barrier height and d the tip- 
surface distance, hence the electric field (~ 1/d) can be 
conveniently tuned through /. By spanning / over 50 
pA— 6 /zA and recording concomitantly the surface state 
spectrum, we establish on an experimental basis the exis- 
tence of a downward Stark shift in .Eo(STS). We also pro- 
vide a one-dimensional model which describes the Stark 
contribution. We extract the Stark-shift-free energy Eq 
of Ag(lll) and find it to agree with the PES value (where 
no electric field is present). 

The Ag(lll) surface was cleaned by Ar + sput- 
ter/anneal cycles, and the measurements were performed 
in a custom-built ultrahigh vacuum STM operating at a 
temperature of 4.6 K, using electrochcmically etched W 
tips, further treated in-situ by indentations into the sur- 
face. Spectroscopy of the differential conductance dl/dV 
versus the sample voltage V was performed by open- 
ing the feedback loop at V = 100 mV in the center of 
impurity- and step-free regions of the surface (> 400 
nm 2 ). A lock- in detection amplifier was employed to 
record the dl/dV (AC voltage modulation was 1 mV rms 
and ~ 1 — 10 kHz, and a single spectrum was acquired in 
~ 3 — 10 s). The dependency of the dl/dV spectrum on 
d is probed through Eq. ^ by varying I. 

Figure ^ summarizes the essential experimental find- 
ings. The sharp step- like rise of the conductance shown 
in Fig. for 7=1 nA allows to determine the low edge 
of the energy band of the Ag(lll) surface state — at 
Eq = —66.7(5) meV in Fig. in agreement with pre- 
vious STS studies 0, but lower than £ (PES)= -63(1) 
meV [7j. In contrast, from the width A of the onset 
we extract T(STS)= 6.4(5) meV 0, in agreement with 
r(PES)= 6.0(5). When increasing I, a sizeable shift of 
Eq occurs. In typical dl/dV spectra (Fig.Qb) the surface 
state onset shifts downward in energy with no apprecia- 
ble broadening as I increases from a pA-range to a fiA- 
range. At / = 6 /iA, the highest current where we could 
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FIG. 2: Eo versus I. The arrow indicates the PES value of 
E 0. 



FIG. 1: (a) Typical dl/dV spectrum taken on Ag(lll) (T = 
4.6 K). The feedback loop was opened at I = 1 nA and V = 
100 mV. A geometrical analysis is used to determine Eo and 
the width A of the surface state 0, 13 • (b) Downward shift of 
Eo when I is increased. All spectra are averages of at least 
5 single spectra from varying sample locations and tips; the 
spectra were renormalized to match the conductance of the 
50 pA spectrum. 



perform reliable spectroscopy, Eq has shifted by —20% 
with respect to J = 50 pA. The dependence of E on /, 
i.e. on the tip-surface distance, is the major finding of 
this Letter. 

Figure [21 presents a quantitative evaluation of Eq for 
all tunneling currents investigated. Two distinct regimes 
are found for the band edge shift — also directly visible 
on the spectra of Fig. ^p: typically 50 pA< I < 1 fiA 
over which Eq decreases by —4 meV, and 2 /iA< I < 6 
/iA where the shift of Eq is more pronounced (—10 meV). 

Care was taken to systematically survey the tip and 
the surface status during spectrum acquisition. We pro- 
ceeded as follows: (i) for a given tip preparation a I = 1 
nA spectrum is acquired and Eq is evaluated following 
the geometrical analysis detailed in Fig.[TJi, (ii) I is then 
set to the value of interest and a spectrum is recorded, 
(iii) the tip and the surface status are checked through 
a final 1 nA spectrum. To minimize errors in the eval- 
uation of the shift, each spectrum acquired in step (ii) 
was fitted to the 1 nA spectrum of step (i) in order to 
determine Eq. For all the data reported in this Letter, no 
change was discernable in the spectra acquired in step (i) 
and in step (iii), meaning that neither a tip modification 
nor a tip-induced damage of the surface occurred when 
acquiring spectra in the 50 pA — 6 [iA range. This con- 
clusion also agrees with the absence of any modification 
in the surface images acquired before step (i) and after 
step (iii). 

Before going any further, a remark concerning the ex- 



perimental setup is necessary. Given the unusually wide 
range of the tunneling resistances R employed in this 
study (from 20 kfl to 2 Gf2), the possibility of a voltage 
drop at the input impedance Ri n of the current pream- 
plifier has to be considered. It would decrease the bias 
voltage at the tunneling junction by (1 + Rin/R) -1 thus 
producing a sizeable i?-dependent, i.e. /-dependent, shift 
of the spectrum when R ~ Ri n . With the variable gain 
current preamplifier employed this does not occur, since 
Rin/R ~ 0.3% at most (this condition is met for both 
DC and AC signals at the modulation frequencies used). 
To highlight this point, Fig. [3^ displays two spectra ac- 
quired at 10 nA but with gains of 10 7 and 10 8 V/A for 
which R m /R = 0.02% and R in /R = 0.1% respectively. 
As expected, their difference is a flat line over the entire 
voltage range (Fig. I3t>), even in the onset region where 
a shift artifact would produce a spike. This is also true 



FIG. 3: (a) dl /dV spectra acquired at 10 nA, 100 mV, with a 
10 7 and a 10 8 V/A gain of the current preamplifier, (b) Differ- 
ences between pairs of spectra acquired at the same currents 
(1 nA, 10 nA, 100 nA and 1 fiA) but with different gains. 
Gains range from 10 5 to 10 9 V/A with R in = 50 fi, 60 ft, 150 
ft, 1 kfl and 10 kfl, respectively. Spectra are shifted vertically 
for clarity. 
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for pairs of spectra acquired at 1 nA, 100 nA and 1/zA at 
different gains. In conclusion, no relevant voltage drop is 
present in our setup at the tunneling currents of interest. 

In order to express the dependency of E on the tip- 
surface distance, we measured the variation of I while ap- 
proaching, with the feedback loop open, the tip towards 
the surface by « 6 A. The recorded currents which cover 
20 pA — 6 [iA follow an exponential behavior up to « 2 
/LtA, from which we extract through Eq. ^ an apparent 
barrier height of <j> = 4.0(2) eV, typical for noble metals 
[IH . In contrast, for / > 2 ^A, I increases faster than 
expected. It is well known from break junction experi- 
ments 0| ! that this signals that the junction is no longer 
in a tunneling regime, rather in a contact regime where, 
as we discuss below, important modifications of the tip 
and the surface morphologies occur. 

Our results are summarized in Fig. where Eq is now 
plotted versus the tip-surface displacement. As shown, 
the negative Stark shift of Eq can be tuned by approach- 
ing or retracting the tip from the surface. To gain further 
insight, we propose a model which describes the exper- 
imental data and enables extrapolation to the zero-field 
properties of the surface state. We first detail some as- 
pects of the model and then discuss the calculation in the 
light of our data. 

The Ag(lll) surface state electrons are modeled with 
the one-dimensional potential proposed by Chulkov et al. 
[l3| . This potential is periodic in the bulk (z < 0), has a 
potential well just outside the surface (0 < z < Z\), then 
decays exponentially towards the vacuum (z\ < z < z- lm ) 
to finally cross over to a long-range image potential 
(z > z; m ). To take into account the presence of the tip 
at z > Ztip, we added to the potential the linear contri- 
bution of the voltage bias V between tip and surface, as 
well as the difference between the work functions of the 
tip (</>tip) and the surface (</> S amp) to include the contact 
potential. Furthermore, we modified the shape of the 
image potential to account for multiple images in the tip 
and the surface. This yields 

V(z) = 2Vm cos gz, z < 

V(z) = V 2 q + V 2 cos pz, < z < zi 

V{z) = V lin (z) + V 3 e-^ z -^\ Zl <z< z im 

V{z) = Vn n {z) - V im (z), z im < z < z tip 

where 

Vn n (z) = E F + s(eV + 4>tip) + (1 - s)0 samp 

where Ep is the surface Fermi energy, s = (z — zi)/(z t ip — 
zi), and 

Vo-KHqKZ-^ ~ Zim) 

with ( — (z — Zi m )/(z^ — zim) and ^ is the digamma 
function. The parameters Vin, g, V20, V2, P, and 



FIG. 4: Eo versus tip-surface distance. Since experimentally 
we only probe a displacement of the tip towards the surface, 
the calculated curves were shifted horizontally to match the 
data (the origin is arbitrarily fixed at 50 pA). Lines show 
calculated results for tip work functions 0ti P : 4.56 eV (solid), 
4.36 eV (dot), 5.2 eV (short dash), and in absence of the tip 
(dash). Also shown, the boundaries to the 2 — 6 fiA region. 



z\ = 57r/(4/3) describing the potential in the absence of 
the tip are fixed to their corresponding values of Rcf. 
The remaining parameters V3, a, A, and z lai and z\^ that 
account for the tip contribution, were fixed by requiring 
the potential and its derivative to be continuous every- 
where, except at z = z t - lp for the derivative. The energy 
of the surface state at a given tip-surface distance was 
calculated by searching for the corresponding maximum 
in the transmission probability for electrons tunneling 
from the tip to the surface (a small imaginary part was 
introduced in the potential between < z < z\ to mimic 
the inelastic scattering). 

The calculation performed with the Ag(lll) work 
function </> sa mp = ^tip = 4.56 eV for both tip and sur- 
face (solid line in Fig. 0} reproduces the data up to / « 2 
fj,A , where a strong discrepancy is observed. We first 
focus on the J < 2 /j,A region were the model works. 
To probe the influence of the contact potential, we per- 
formed calculations with different values for (j> t i p : for a 
pure W tip (^ t i P = 5.20 eV, short dashed line on Fig.^J, 
and for a Ag coated tip with a high density of surface de- 
fects ((fitip = 4.36 eV, dotted line). The contact potential 
yields, in the former Stark shift opposite to the 

one of the negative bias, whereas the latter increases the 
downward shift. As shown, neither of them reproducse 
the experimental behavior as well as the calculation with 
0ti P = 4.56 eV does. Although we cannot fully exclude a 
contribution from the contact potential, from these con- 
siderations we conclude that, for the tips used, it must 
bear a negligible contribution to the Stark effect. We 
note also that since the data was obtained with a variety 
of tips, the highly reproducible shift observed indicates 
that our in-situ tip preparation must lead to similar tip 
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apices. In conclusion, on the basis of our model, at the 
largest tip-surface distances shown, the Stark effect is 
mainly produced by the bias and by the lowering of the 
image potential compared with the case of an isolated 
surface, with a 2 to 1 ratio between the two. However, as 
tip-surface distance decreases the lowering of the image 
potential is increasingly stronger and eventually governs 
the Stark effect. 

Now we discuss the asymptotic behavior of the Stark 
effect. As the tip is retracted, the Stark contribution 
decreases towards the non-perturbed Eq binding energy 
which is calculated by suppressing the presence of the 
tip in our model. It yields a constant value (dashed 
line on Fig. with tip-surface distance of Eq = —64 
meV, in agreement with the Stark-free PES value for 
Ag(lll). Figure 01 shows that the Stark effect cannot be 
suppressed, as the tip would have to be retracted well be- 
yond the lowest tunneling currents experimentally acces- 
sible, which are typically in the pA-range. For Ag(lll) 
the Stark shift results then, at the best, in a ~ 4 — 5% 
error in the evaluation of Eq. 

We now turn to the / > 2 \xA region where our model 
calculation no longer reproduces the experimental data. 
The comparison between the two indicates that the ac- 
tual electric field in the junction is stronger than the one 
predicted in our model calculation. Since in the model 
the tip and surface morphologies are assumed to be con- 
stant at all tunneling currents, we hint that in the I > 2 
/iA range, i.e. R < 50 kO, the tip and the surface deform 
to yield a stronger electric field. However, this deforma- 
tion is reversible, since the spectra acquired at / = 1 nA 
prior and after ramping the current to I > 2 fiA (three- 
step procedure described earlier) did not yield any sub- 
stantial differences. This picture agrees with calculations 
performed for Au(100) which indicate that the tip 
and the surface undergo an elastic deformation, i.e. they 
stretch towards each other, because of attractive adhe- 
sive forces acting at R < 100 kfi. In particular, this is 
predicted to result in a deviation from the exponential 
behavior of Eq. ^ which we do indeed observe in the I 
versus displacement measurements at I > 2 fiA. 

To summarize, in STS a downward shift of Eq of the 
Ag(lll) Shockley surface state occurs when decreasing 
tip-surface distance. A model calculation explains this 
observation in terms of a Stark shift produced by the bias 
and by the lowering of the image potential compared to 
an isolated surface. We extract a non-perturbed value of 
—64 meV for the binding energy of the Ag(lll) surface 
state. In light of the agreement with the PES measure- 
ments for Ag(lll), we conclude that a possible positive 
shift of Eq in the PES data of Ref.tD due to surface defects 
is negligible. Finally, for tunneling resistances R < 50 
kO, an enhanced shift is observed, which we assign to an 
elastic deformation of the tip-surface morphology. These 
results, supported by our recent observation of a Stark 
effect for other noble metal surface states and for the 



Ag/Ag(lll) adatom [3, suggest that STS data require 
a thorough quantification of the Stark effect when striv- 
ing for high energy resolution, especially for states whose 
wavefunctions have large decay lengths into vacuum. 
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